Anisotropic electron-phonon interaction in the cuprates 
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We explore manifestations of electron-phonon coupling on the electron spectral function for two 
phonon modes in the cuprates exhibiting strong renormalizations with temperature and doping. 
Applying simple symmetry considerations and kinematic constraints, we find that the out-of-plane, 
out-of-phase O buckling mode (-Big) involves small momentum transfers and couples strongly to elec- 
tronic states near the anti-node while the in-plane Cu-O breathing modes involve large momentum 
transfers and couples strongly to nodal electronic states. Band renormalization effects are found to 
be strongest in the superconducting state near the anti-node, in full agreement with angle- resolved 
photoemission spectroscopy (ARPES) data. 

PACS numbers: 



The electron-phonon interaction in metals is generally 
considered to be isotropic. Recent developments in the 
cuprates however raises important questions about this 
assumption. The renewed interest stems from the obser- 
vation of bosonic band renormalization effects in a wide 
range of p-type materials [1]. This effect is found to be 
dramatically enhanced below T c near the (ir, 0) regions in 
the Brillouin zone (BZ)[2]. Such a momentum and tem- 
perature dependence appears to support the idea that the 
bosonic renormalizations in p-type materials is due to a 
coupling of the electrons to the 41 meV spin resonance 
mode that appears below T c [3, 4]. 

On the other hand, it has been argued that the bosonic 
renormalization effect may be reinterpreted as a result 
of a coupling to the half-breathing in-plane Cu-0 bond- 
stretching phonon for nodal directions [5] and the out- 
of-plane out-of-phase O buckling B\ g phonon for anti- 
nodal directions [6]. This re-interpretation has the clear 
advantage over the spin resonance in that it naturally 
explains why the band renormalization effect is observed 
under three circumstances: 1) in materials where no spin 
mode has been detected, 2) in the normal state, and 3)in 
the deeply overdoped region where the spin mode is nei- 
ther expected nor observed. However this interpretation 
also requires that the electron-phonon interaction to be 
highly anisotropic and its impact on the electrons to be 
strongly enhanced in the superconducting state. This is 
something one does not expect a priori. 

This paper investigates the anisotropy of the electron- 
phonon interaction of two known phonons - the half- 
breathing mode and the B\ g buckling mode - formu- 
lated within the same framework as that used to explain 
phonon lineshapc changes observed via Raman and neu- 
tron measurements [4, 7, 8, 9]. We find that the electron- 



phonon interactions to be very anisotropic as a conse- 
quence of the following four properties of the cuprates: 
1) the symmetry of the phonon polarizations and elec- 
tronic orbitals involved leading to highly anisotropic 
electron-phonon coupling matrix elements g(k, q); 2) the 
kinematic constraint related to the anisotropy of the 
electronic band structure and the van Hove singularity 
(VHS); 3) the d-wave superconducting gap, and 4) the 
near degeneracy of energy scales of the B\ g phonon, the 
superconducting gap, and the VHS. Taking into account 
these and thermal effects we can naturally explain the 
observed momentum and temperature dependence of the 
data, leading to unified understanding of band renormal- 
izations. 

Our starting point is the tight-binding three-band 
Hamiltonian modified by B\ g buckling vibrations as con- 
sidered in Ref. [10] and in-plane breathing vibrations as 
considered in Ref. [11]: 

H = ^ ec u bt,a b n,<7 

+ >, Ps[t- Qsgdpuf (njK^/n,^ + h.c.) 
+*' Ps,6'(al iS a n ,s',„ + h.c) (1) 



n,5,5> 



with Cu-0 hopping amplitude t, 0-0 amplitude i? , and 
ed, p denoting the Cu and O site energies, respectively. 
Here a a=x _ y , oJ a=x y (b, tf) annihilates, creates an electron 
in the O a (Cu) orbital, respectively. The overlap factors 
Q$,P$, and Ps,s' are given in our notation as Q± x = 
Q± y = ±i;-P±x = -P±y = ±1;-P±x,±j/ = 1 = —P?x,±y 
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For the half-breathing mode we consider only couplings 
9dp — qot (with 1/qo a length scale) arising from modu- 
lation of the covalent hopping amplitude t [11, 12]. For 
the Big mode a coupling to linear order in the atomic 
displacements arises from a local c-axis oriented crystal 
field E z which breaks the mirror plane symmetry of the 
Cu-0 plane[10]. 

The symmetry of the p, d wavefunctions of the Cu- 
O plane is embedded in the amplitudes </> of the tight- 
binding wavefunctions forming the anti-bonding band 
with energy dispersion e(k) arising from Eq. 1. Specifi- 
cally the transformation is determined by the projected 
part of the tight-binding wavefunctions = <p b {k)d^ a 
and a a ,k,a = Q (k)d k , CT , 
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N(k) 



[e 2 (k)-t' 2 (k)], 



(2) 



,(k) 



AT(k) 1 



[e(k)t x , y (k)-t'(k)t y , x (k)], (3) 



with A 2 (k) = [e 2 (k)-t' 2 (k)] 2 + [e(k)t x (k)-t'(k)t y (k)} 2 + 
[e(k)t y (k) - t'(k)^(k)] 2 , t a (k) = 2t sin(fc Q a/2), and 
t'(k) = At' sm{k x a/2) s\n{k y a/2). 

The projected wavefunctions act in conjunction with 
the phonon eigenvectors to determine the full sym- 
metry of the electron-phonon coupling: H e i- p h — 

7f Ek,q,,,„ 3m ( k , q)4, CT c k+q,^ [4,q + «m, -q] ■ Neglecting 
the motion of the heavier Cu atoms compared to O in the 
phonon eigenvectors the form of the respective couplings 
can be compactly written as 



9B lg (k, q) = eE z 



_ 4M M(q)0 Bl 
x {<f) x (k)<j) x (k') cos(q y a/2) - <j> y (k)</> y (k') cos(g x a/2)} 

g br (k,q) = g dp 



(4) 



2M tt 



br 



E 



(5) 



x {0 6 (k')</> a (k) cos(fc>/2) - (j) b {k)<f> a {k') cos(fc Q a/2)} , 

withM(q) = [cos 2 (q x a/2)+cos 2 (q y a/2)}/2 and k' = k-q. 
A simple consequence of symmetry can be observed at 
small momentum transfers q where gs lg {k, q = 0) ~ 
cos(fc x a) — cos(fcj,a), while g br ~ sin(ga) for any k. Yet 
for large momentum transfer q = (w/a,Tv/a), gB lg van- 
ishes for all k while g br has its maximum for k located on 
the nodal points of the Fermi surface. This g-dependcncc 
has been the focus before in context with a d x i_ y i pair- 
ing mechanism [13, 14], but the dependence on fermionic 
wavevector k has usually been overlooked. It is ex- 
actly this fermionic dependence which we feel is crucially 
needed to interpret ARPES data. 

These features of each coupling constant can be seen 
in Fig 1 which plots | g(k, q) | 2 as a function of phonon 
scattered momentum q connecting initial k and final 



FIG. 1: Plots of the electron-phonon coupling | g(k, q) | 2 for 
initial k and scattered k' = k — q states on the Fermi surface 
for the buckling mode (left panels) and breathing mode (right 
panels) for initial fermion k at an anti-nodal (top panels) 
and nodal (bottom panels) point on the Fermi surface, as 
indicated by the arrows. The red/blue color indicates the 
maximum/minimum of the el-ph coupling vertex in the BZ 
for each phonon. 



k' = k q fermion states both on the Fermi surface. 
Here we have adopted a five-parameter fit to the band 
structure used previously for optimally doped Bi-2212[3]. 
For an anti-nodal initial fermion k^^v, the coupling to the 
Bi g phonon is largest for 2kp scattering to an anti-nodal 
final state. The coupling to the breathing mode is much 
weaker overall, and is especially weak for small q and 
mostly connects to fermion final states at larger q. For a 
nodal fermion initial state kjv, the Bi g coupling is sup- 
pressed for small q and has a weak maximum towards 
anti-nodal final states, while the breathing mode attains 
its largest value for scattering to other nodal states, par- 
ticularly to state — kjv. 

In Fig. 2 we show effective couplings Ak (summed over 
all q) for k points in the first quadrant of the BZ in the 
normal state[14] . Here we have used the value of the elec- 
tric field eE z = 1.85eV/ A consistent with Raman scatter- 
ing measurements of the Bi g phonon and buckling of the 
Cu-0 plane in YBaCuO, and have taken qot = 2^/2eE z 
as a representative value for the coupling to the breathing 
mode. These values lead to net couplings similar to that 
obtained from LDA calculations[14]. As a result of the 
anisotropy of the coupling constant, large enhancements 
of the coupling near the anti-node for the Bi g phonon, 
and nodal directions for the breathing phonon mode are 
observed, respectively, which are an order of magnitude 
larger than the average A's. 

We now turn to the spectral function. The dynamic 
nature of the coupling is governed by the energy scales 
of the band structure, phonon frequencies, and super- 
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FIG. 2: Plots of the electron-phonon coupling Ak in the first 
quadrant of the BZ for the buckling mode (right top panel) 
and breathing mode (right bottom panel). The color scale is 
shown on the right for each phonon. The left panel shows 
energy contours for the band structure used [3]. 



conducting gap energy, and are manifest in features of 
the Nambu-Eliashberg electron-phonon self-energy [15]: 
E(k,iw m ) = iu m (l - Z(k,iuj m ))T + x(k, iu m )f 3 + 
3>(k, iuj m )ri with 



£(k,wj m ) = — fM( k 'P- k ).9p(P ? k -P) 

xD° (k-p, iuj m - iiv n )f 3 G°(p,iuj m )T 3 , 



(6) 



with Gr°(p, iu>, 
i a 2 



) = ^KSsixr 1 and E2 w - 

e"(k) + A 2 (k), with fj = rj...3 are Pauli matrices. We take 
A(k) = Ao[cos(fca;a) — cos(k y a)]/2 with Ao = 35meV. 
The bare phonon propagators are taken as (q, = 

\ in -n in +n ■ Here we take dispersionless opti- 
cal phonons hfls lg ,br = 36, 70meV, respectively. 

The spectral function A(k,w) = — ^ImGii(k,iuj — > 
lu + iS) is determined from Eq. 6 with G{k, iuj m ) = 
[T iuj m ZCk, iuj m ) - (e(k) + x(k, iu m ))f 3 - $(k, iw m )fi] _1 
and the coupling constants given by Eq. 4-5. Here we 
consider three different cuts: along the nodal direction 
(k x = k y , denoted by c), and two cuts parallel to the 
BZ face (k y = 0.757r/a, b) and (k y — 0.647r/a, a) where 
bilayer effects are not severe [6], as indicated in Fig. 2. 

Our results for the three cut directions are shown in 
Fig 3 for the normal (T = 110K) and superconduct- 
ing (T = 10K) states. One can only observe very weak 
kink features in the normal state, but a much more pro- 
nounced kink at 70 meV in the superconducting state is 
observed for all the cuts, in excellent agreement with the 
data[6], shown also in Fig. 3. The dramatic tempera- 
ture dependence stems from a substantial change in the 
electronic occupation distribution and the opening of the 
superconducting gap. In the normal state, the phonon 
self energy is a Fermi function at 110K centered at the 
phonon frequency, which results in a thermal broaden- 
ing of 4.4ksT or 50meV, comparable to phonon frequen- 
cies. At 10K in the superconducting state, on the other 
hand, the phonon self energy is sharply defined due to 



the step-function-like Fermi function and a singularity 
in the superconducting density of states. These factors, 
together with the experimental difficulty in detecting a 
weak effect, conspire to leave the false impression that 
the bosonic mode merges below T c , as invoked by the 
spin resonance scenario [3]. 

In the superconducting state, the data and the theory 
exhibit a rich variety of signatures in the electron-phonon 
coupling: cut a exhibits the strongest signatures of cou- 
pling, showing an Einstein like break up into a band that 
follows the phonon dispersion, and one that follows the 
electronic one; in cut b, a weaker signature of coupling 
manifests itself in an s— like renormalization of the bare 
band that traces the real part of the phonon self-energy; 
and in cut c, where the coupling to the B\ g phonon is 
weakest, the electron-phonon coupling manifests itself as 
a change in the velocity of the band. The agreement 
is very good, even for the large change in the electron- 
phonon coupling seen between cuts a and b, separated by 
only l/10th of the BZ. 

We now discuss the reason for the strong anisotropy as 
a result of a concurrence of symmetries, band structure, 
d— wave energy gap, and energy scales in the electron- 
phonon problem. In the normal state, weak kinks are ob- 
served at both phonon frequencies, but more pronounced 
at 70 meV for cut c and at 36 meV for cut a, as ex- 
pected from the anisotropic couplings shown in Fig. 1. 
However this is also a consequence of the energy scales 
of the phonon modes and the bottom of electronic band 
along the cut direction as shown in Fig. 2. Further to- 
wards the anti-node, another interesting effect occurs. As 
one considers cuts closer to the BZ axis, the bottom of 
the band along the cut rises in energy and the breath- 
ing phonon mode lies below the bottom of the band for 
(k x = 0, k y > 0.827r/a) and so the kink feature for this 
mode disappears. This is the usual effect of a Fano re- 
distribution of spectral weight when a discrete excitation 
lies either within or outside of the band continuum[16]. 
This thus can be used to open " windows" to examine the 
coupling of electrons to discrete modes in general. Taken 
together, this naturally explains why the normal state 
nodal kink is near 70 meV[5] while the anti-nodal kink is 
near 40 meV[2, 6]. 

In the superconducting state the bare band renormal- 
izes downward in energy due to the opening of the gap 
and gives the strongest effect for the anti-nodal region 
where the gap is of the order of the saddle point en- 
ergy. This also gives the appearance of an increase of 
bandwidth. The B\ g coupling becomes dramatically en- 
hanced by the opening of the gap at a frequency close to 
the frequency of the phonon itself, and leads to the dra- 
matic renormalization of the phonon observed in Raman 
and neutron measurements [7]. Yet the breathing cou- 
pling does not vary dramatically as the large gap region 
is not weighted as heavily by the coupling constant. As a 
consequence, the B\ g mode dominates and clear kinks in 
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FIG. 3: Image plots of the calculated spectral functions in the normal (al,bl,cl) and superconducting (a2,b2,c2) 
states compared to the spectral functions in the normal (a3,b3,c3) and superconducting (a4,b4,c4) states measured in 
Bi2Sr2Cao.92Yo.o8Cu208+5 (Bi-2212)[6] for momentum cuts a,b,c shown in the right-most panel and in Figure 2. The same 
color scale is used for the normal/superconducting pairs within each cut, but the scaling for the data and the calculation are 
separate. The red markers indicate 70meV in the superconducting state. 



the data are revealed at around 70 meV - the B\ g energy 
plus the gap. Indeed a much weaker kink at 105 meV 
from the breathing phonon shows at higher energies for 
the nodal cut but this kink becomes weaker away from 
nodal directions as the coupling and the band bottom 
along the cut moves below the phonon frequency so that 
the "energy window" closes. 

Our work is formulated in the context of band struc- 
ture and one phonon process. Fig. 3 indicates that 
such an approach catches the key features of the band 
renormalization effects due to optical phonons in the 35- 
70 meV range, something unexpected a priori given the 
strong correlation effects in these materials. Additional 
broadening in the data of Fig. 3 is likely a manifestation 
of the Coulomb correlation over a broad energy range. A 
next step is to incorporate the correlation effect, vertex 
correction, and possible non-linear phonon effects, espe- 
cially as we push towards underdoped regime. These are 
challenges ahead. On the other hand, the current ap- 
proach allows us to study the doping trend in the over- 
doped regime. 

In summary, contrary to usual opinion of the role of 
anisotropy in electron-phonon coupling, the interplay of 
specific coupling mechanisms of electrons to the buck- 
ling and breathing phonons give a natural interpretation 
to the bosonic renormalization effects seen in ARPES in 
both the normal and superconducting states, and pro- 
vides a framework to understand renormalizations as a 
function of doping. 
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